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Mean-Field Description of Fusion Barriers with Skyrme's Interaction 
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Fusion barriers are determined in the framework of the Skyrme energy-density functional together 
with the semi-classical approach known as the Extended Thomas-Fermi method. The barriers 
obtained in this way with the Skyrme interaction SkM* turn out to be close to those generated 
by phenomenological models like those using the proximity potentials. It is also shown that the 
location and the structure of the fusion barrier in the vicinity of its maximum and beyond can be 
quite accurately described by a simple analytical form depending only on the masses and the relative 
isospin of target and projectile nucleus. 
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I. INTRODUCTION 

The knowledge of the collective potentials between two 
colliding ions is absolutely crucial for the synthesis of new 
isotopes. This problem has been the subject of a very ac- 
tive research over the last decade and remains one of the 
most intensively studied subjects in low-energy nuclear 
physics in particular in the perspective of the synthesis of 
super-heavy elements as well as of exotic nuclei far away 
from the /^-stability line. It has in particular been shown 
that models based on a macroscopic approach such as 
the Liquid Drop Model (LDM) or of semi-classical type 
like the Extended Thomas-Fermi (ETF) method together 
with the Skyrme energy-density functional are able to 
reproduces quite accurately experimental data on fission 
and fusion barriers. 

The research concerning the interaction potential be- 
tween colliding ions goes back to the work of R. Bass (for 
a review see ref. 0). Since the idea of a proximity po- 
tential due to W. J. Swiatecki and coworkers in the late 
1970's 0, many improvements have been proposed to 
make this phenomenological approach more realistic and 
general, in particular by taking into account the local 
curvature of the surfaces of target and projectile 0, Q . 
One of the main challenges was, as already mentioned, 
to give a reliable guideline for the formation and stability 
of super-heavy elements. 

The parameters of the proximity function are usually 
fitted to experimentally known fusion barriers heights. 
The basic idea of all proximity models is to determine 
the potential between the two colliding nuclei as func- 
tion of the minimal distance s of their surfaces and their 
so-called reduced radius defined in the case of spherical 
nuclei as R = R\ R%/ '(R1+R2), where the indices 1 and 
2 refer to the target and projectile nucleus respectively. 

Experimental data confirm the existence of a pocket 
in the shape of the nuclear part of the potential. This 
feature plays an essential role for the formation proba- 
bility but also for the stability of the compound nucleus 
in the fusion experiments. Indeed, if the pocket is deep 
and wide, several quasi-bound states might exist and the 
probability of forming that compound nucleus is large. 



If, on the other hand, this minimum is shallow and nar- 
row no such states will exist. It would therefore be of 
great help to have at hand a simple, yet sufficiently accu- 
rate phenomenological expression for the determination 
of fusion barriers in order to evaluate the most favorable 
reaction for the formation of a given compound nucleus. 

The aim of our present investigations is to show that 
the shape of the proximity potential can be, indeed, 
quite accurately determined through self-consistent semi- 
classical calculations. We proceed in a way similar to that 
suggested in Ref. [j| and recently used in Ref. 0. Af- 
ter introducing the Skyrme Hartree-Fock energy-density 
functional together with the semi-classical ETF approx- 
imation in Section II we explain in Section III how this 
approach can be used to calculate fusion barriers. In 
Section IV the fusion barriers obtained in this way for 
a sample of 269 different reactions (different combina- 
tions of target and projectile) leading a priori to different 
isotopes of the super-heavy elements with even charge 
numbers in the range Z = 108 — 114 are presented and 
compared to the barriers obtained in other theoretical 
models. Based on these fusion-barrier calculations we 
present in Section V a simple analytical expression for 
this potential and show that it is able to reproduce very 
accurately the fusion barrier heights obtained in the ETF 
approach. Conclusions are finally drawn in Section VI. 



II. THE SKYRME ETF APPROACH TO 
NUCLEAR STRUCTURE 

For the calculation of the interaction potential be- 
tween the two nuclei we need an energy-density func- 
tional which has proven its capacity to correctly describe 
nuclear ground-state properties. As we are interested 
in the potential energy resulting from the interaction of 
the tails of the density distributions of target and pro- 
jectile, a very precise description of the nuclear surface 
seems absolutely crucial. The Skyrme effective interac- 
tion has been shown to have this pro pert y 0, H, @ an d 
has been, together with Gogny's force |10|. the most suc- 
cessful effective nucleon-nucleon interaction over the last 
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three decades. Its energy density is given by 
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where the constants Bi are combinations of the usual 
Skyrme- force parameters tj and Xj . ^ r q and J q des- 
ignate the kinetic energy density and spin-orbit density 
of charge state q, q — {n,p} and non indexed quantities 
correspond to the sum of neutron and proton densities 
such as p = p n + pp. 

Since our aim is to evaluate fusion barriers, a quan- 
tity which is essentially determined by the radii and tails 
of the density distributions of target and projectile, we 
have chosen the semi-classical approximation known as 
the Extended Thomas-Fermi (ETF) method 01 to de- 
termine in a self-consistent way the structure of projectile 
and target nuclei. 

The ETF approach allows to express T q and J q as func- 
tions of the local density p q and its derivatives. This 
functional expression has the form of an expansion with 
Planck's action quantum h as order parameter. For prob- 
lems which don't break time-reversal symmetry (as would 
e.g. be the case for rotating nuclei) only even powers of H 
appear so that the semi-classical expansion of the kinetic 
energy density will take the form 
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At lowest order in the ETF expansion one has the well- 
known Thomas-Fermi kinetic energy density 
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whereas the second-order contribution is given by 
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Here the effective- mass form factor f q = m/m*(r), is 
given by 



f q (r) 



m 



2m f 5£{r) 



m*(r) h 2 \5T q (r) 
2m 



(5) 



1 



h 2 



[B 3 p(r) + B 4Pq (f)] 



and the spin-orbit potential as 
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In the ETF approach one also obtains a semi-classical 
expansion for the spin-orbit density. Since the spin is 
a pure quantal phenomenon which does not have any 
classical analog, the semi-classical expansion of J starts 
at order h 2 and the leading term is given by 
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This now, however, means that one is able to express 
the total energy density £ as a unique functional of the 
proton and neutron density alone and to perform a vari- 
ational calculation of the ground-state properties of a 
nucleus, where the variational quantities are p n (r) and 
p n (r ) . That such a unique functional must exist for a sys- 
tem of interacting Fermions has been shown by Hohcn- 
berg and Kohn 12]. In the general case this functional 
is, however, perfectly unknown. It is the achievement of 
the ETF approach to derive this functional in the semi- 
classical limit in a systematic way. ETF variational cal- 
culations have indeed proven, when the semi-classical ex- 
pansion is carried through to order H 4 |lllll3| . extremely 
successful to describe average (liquid-drop type) nuclear 
properties (see Ref. QU| for a review). It has in particular 
been shown 0, IliL Il5| that modified Fermi functions of 
the form 
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are in that case a very good approximation to the full 
variational solution. A systematic analysis of the varia- 
(3) tion of the density parameters p^\ a q and j q (the pa- 
rameter R q is then determined through particle-number 
conservation) which are the variational parameters of the 
ETF approach, with the mass number A and the isospin 
parameter / can be found in Ref. |16( . 

We have carried out these kind of 4 th order ETF calcu- 
lations, to determine the structure of target and projec- 
tile nucleus. The densities obtained in this way are then 
used to determine the fusion potential between the two 
nuclei as will be explained in the following. 



III. THE ETF FUSION POTENTIAL 

For evaluating the fusion potential we will use the sud- 
den approximation, i.e. keeping the densities of the collid- 
ing ions fixed and neglecting all possible rearrangement 
effects. At first sight, this looks like a quite crude ap- 
proximation. For a fusion reaction the beam energy per 
nucleon is, indeed, quite small as compared to the Fermi 
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energy and an adiabatic treatment would be more appro- 
priate. One has, however, to keep in mind that the essen- 
tial characteristics of the fusion barrier, as its structure 
around the above mentioned pocket, are determined at 
distances larger and close to the touching configuration. 
Defining the minimal distance s between the equivalent 
sharp surfaces of the liquid drops identifying the two nu- 
clei we see that the essential features of the fusion barrier 
or cross-sections are determined at the distances (s > 0), 
where the sudden approximation seems to be a very rea- 
sonable approach. We can then evaluate the nuclear part 
of the interaction potential V nuc of two spherical colliding 
ions as function of the distance d between theirs centers 
of mass (see Fig. 1) as : 
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FIG. 1: Definition of target and projectile coordinates. 



V nuc (d) = / \£ nuc {pV(r) + pW(f- d)} 



(9) 
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where £ nuc is the nuclear part of the energy-density, Eq. 
||TJ, and pW(r') denotes the density distribution of nu- 
cleus i alone ((1) labeling e.g. the target and (2) the 
projectile). 

If the distance d between the two nuclei is much larger 
than the sum of the half density radii (d ^> R 1 +R 2 ) then 
the nuclear part of the interaction should vanish because 
of the short range of the nuclear forces: 
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Thus, the nuclear part of the interaction potential i|l(J|) 
between the nuclei is non negligible only when the dis- 
tance s between the nuclear surfaces of both nuclei is 
smaller then a few fm. The Coulomb interaction between 
the two ions is given as usual in the form 



V Co ui(d) = 



\ri~r 2 \ 



d 3 ri d 3 
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where we simply assume that the charge density is given 
in terms of the proton density by p^ w e pp ' where the 
latter can be approximated in the ETF approach by Eq. 

©■ 



IV. FUSION BARRIERS FOR SUPER-HEAVY 
ELEMENTS 

The fusion barrier appears as the result of the compe- 
tition between the long-range repulsive Coulomb interac- 
tion Vcoui and the short-range attractive nuclear forces 
V nuc as is illustrated on Fig. 2, where the total fusion bar- 
rier is plotted as function of the distance d (see Fig. 1) of 
the centers of mass of projectile and target nucleus as well 
as of the distance s between their half-density surfaces. 
One notices that for distances s > 3 fm the attractive 
nuclear part of the fusion potential almost vanishes and 
the barrier is just determined by the repulsive Coulomb 
potential. Taking the diffuseness of the charge-density 
distribution into account reduces the height of the fusion 
barrier only by about 1 MeV as compared to the barrier 
obtained with point charge distributions of target and 
projectile. 




FIG. 2: Shapes of total (nuclear + Coulomb) fusion barriers 
V (full line) as function of the center-of-mass distance d (top) 
and the distance between the equivalent sharp surfaces s (bot- 
tom) for the reaction 50 Ti + 238 U leading to the super-heavy 
element Z = 1U. Also shown are the exact Coulomb barrier 
(thin line) and the one corresponding to two point charges 
(dashed line). 

It is evident that when the overlap of the densities of 
the colliding nuclei becomes large, as e.g. when s < 
(dotted line in Fig. 2) the sudden approximation used 
here becomes more and more questionable. Our aim, 
however, is not to give a precise description of the entire 
fusion potential, i.e. also for large negative s values, but 
rather to obtain some reasonable estimates of the height 
Eb of the fusion barrier. Fig. 3 shows Eb obtained in 
the sudden approximation and neglecting deformation ef- 
fects of the colliding ions for 269 different reactions (with 
/3-stable targets ranging from 168 Er to 238 U) leading to 
different isotopes of the super-heavy elements with even 
charge numbers between Z = 108 and Z = 114. The 
barriers shown are the ones obtained with the Skyrme 
forces SkM* Q (barriers corresponding to the Skyrme 
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FIG. 3: Height of fusion barriers corresponding to the different reactions leading to different isotopes of the S.H. elements 
with Z=108, 110, 112, 114, evaluated in our ETF approach (thick line) and in the proximity model of Ref. [T^ (thin lines), as 
function of the projectile neutron number. 



SLy4 interaction Q are identical to within 1 MeV). The 
thin lines illustrate the heights of the fusion barriers ob- 
tained in the Myers-Swiatecki proximity model [l7| . Bar- 
rier heights evaluated using the prescription of Ref. 
are very close to the former differing by not more than 
2 MeV for all the reactions considered here. 



V. A SIMPLE ANALYTICAL FUSION 
POTENTIAL 

The nuclear part of the potential between two colliding 
ions can be quite well reproduced by the approximate 
form 



V nuc (d)*V nuc (d)=V e-^ d - d ^ 



(12) 



where d is the distance between the centers of mass of 
target and projectile and do the distance for which the 
nuclear potential has its minimum, do and the depth Vq 
of the potential are a direct result of our semi-classical 
calculation. These two quantities are parametrized in the 
following way as functions of the masses A\ , A 2 and re- 
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where the values of the parameters entering into these 
expressions can be adjusted to the ETF nuclear poten- 
tials. 

To obtain the potential that the system of two colliding 
nuclei experiences we have to add the Coulomb potential, 
Eq. 1)11(1 . to the nuclear part. For distances large com- 
pared to the sum of the equivalent-sharp-surface radii, 
this Coulomb potential is very well approximated by the 
Coulomb potential between two point charges, but as 
soon as the densities of the two nuclei interpenetrate 
this approximation will no longer give good results and 
it turns out that the difference 



&V Coul (d) = V Coul {d) 



(15) 



5 



will be, in the case of the super-heavy elements consid- 
ered here, of the order of 3 MeV already in the vicinity 
of the touching configuration (s = 0). 

Instead of calculating the diffuse-surface Coulomb en- 
ergy explicitly for each reaction under consideration, the 
idea is to develop an easy-to-use expression to approxi- 
mate in a reliable way the ETF fusion potential between 
two nuclei developed above, including the full Coulomb 
potential. That is why we propose to write this potential 
in the form 



V(d) = V nuc (d) + 



(16) 



where the term V nuc , which now contains the Coulomb- 
energy difference AV Coul , Eq. (fl5|l . in addition to the 
nuclear potential (|1U|) . is approximated by Eq. I|12|) with 
the parameters do an d Vb being given by Eqs. i|13|) and 
(|14(l . It is clear in this context that these new parameters 
can turn out to take values that could be quite different 
as compared to those describing the nuclear part alone. 

To determine these quantities we perform a simultane- 
ous fit of the 5 free parameters to the fusion potentials of 
the 269 reactions obtained in the above described ETF 
approach (see Fig. 3). We proceed in the following way: 
What we are interested in is primarily the precise de- 
scription of the location d max and the height V max of 
the total fusion barrier, and we therefore attach a max- 
imum weight to their best possible reproduction. When 
going to smaller distances, i.e. towards negative values 
of s, the Coulomb-energy difference AV Coul , Eq. (jT5|l . 
grows larger, but also the nuclear potential obtained in 
our sudden approximation becomes less and less reliable. 
That is why we attach a decreasing importance to the re- 
production of the fusion barriers for smaller and smaller 
distances. For increasing distances d > d maxl on the 
other hand, both the nuclear potential V nuc as well as 
the Coulomb correction AV Coul go rapidly to zero and 
our approximation form (|16fl becomes better and better. 
We therefore attach again less weight to the precise re- 
production of the values with large positive values of s. 
We have thus imagined a least-square-fit procedure with 
a normalized weight function of Gaussian form and with 
a width of 1 fm centered at d — d max . The values of the 
parameters obtained in this way are listed in Table 1. 



r [fm] 


b [fm] 


v [MeV] 


K 


a [fm 2 ] 


1.183 


-2.400 


-46.07 


-0.4734 


0.173 



Tab. 1: Values of parameters entering the fit of the 
fusion potential, Eq. i|16|) 

The final criterion of the accuracy of our approach con- 
sists in its ability to reproduce the height and shape of 
the fusion barriers for all 269 reactions obtained within 
the ETF approach with the SkM* Skyrme interaction. 
The r.m.s. deviation of the exact ETF fusion potential 
and its approximation by Eq. (|16f) is only of 0.27 MeV at 



the top of the barrier and of 0.37 MeV for the touching 
configuration (s = 0). 

Such a fusion barrier is shown on Fig. 4 for the re- 



actions 8 Ce 



-Th. As one can see, the ETF barrier 



is almost perfectly reproduced by our analytical expres- 
sion, Eq. I|16fl with (|12|) . The barrier height is also almost 
the same as that obtained in the proximity model of Ref. 
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FIG. 4: Fusion barrier for the reaction 48 Ca + 232 Th obtained 
within the ETF approach (full line), its approximate analyt- 
ical form, Eq. (dashed line) and the proximity potential 
of Ref. E3 (dotted line). 

It is obvious from Fig. 3 that the agreement between 
our ETF approach and the proximity model is not always 
going to be that close. If one considers e.g. the reaction 
110 Cd + 154 Sm one already concludes from Fig. 3 that 
the barrier heights of the two approaches are going to 
be different by some 7 MeV (or 2%). One also notices 
that the minimum obtained in our ETF approach is quite 
shallow, whereas the proximity potential does not make 
any predictions about the region s < 0. One could now 
speculate about the validity of our approach close to and 
inside the touching configuration. The enhancement of 
the height of the fusion barriers related to the sudden 
approximation is obviously related to the concept of the 
diabatic fusion barriers advocated by W. Norenberg and 
co-workers jl9j . 

The coincidence of the ETF fusion barrier heights with 
those obtained in the proximity approach for very asym- 
metric reactions (Ai » A 2 ) seems to be related to the 
value s max of the distance s of the equivalent sharp sur- 
faces of the two colliding nuclei at the top of fusion bar- 
rier, as this is demonstrated on Fig. 5. Comparing Figs. 
3 and 5 one can, indeed, observe that the enhancement 
of the ETF barrier height relative to prediction of the 
proximity model is larger as the distance of the sharp 
surfaces of the two nuclei at the top of the barrier gets 
smaller which is the case for the reactions between nearly 
symmetric nuclear systems. 

As we have just considered fusion reactions leading 



6 



s max = (2.261 - 0.000537 Z^Z 2 ) fm 




500 1000 1500 2000 2500 3000 3500 
Z 1 Z 2 

FIG. 5: Location s max of the top of the fusion barrier as 
function of the product Z 1 Z 2 of the charge numbers of target 
and projectile. 

to the super-heavy elements with even values of Z be- 
tween Z — 108 — 114 one might doubt about the utility 
of such an analysis for the description of fusion barriers 
in other regions of the atomic chart. We have therefore 
applied our phenomenological expression to the reactions 

160,162,166 Dy + 46-50 Ti ag weU ag tQ the very light syg _ 

tern 48 Ca + 48 Ca. To our great surprise the predictions 
based uniquely on our knowledge of the ETF barriers of 
the super-heavy system gave astonishingly good results 
for these lighter systems with a deviation between ETF 
results and those of the simple analytical expression, Eq. 
H16| l. of less than 0.7 MeV for the former and of about 
1.3 MeV for the latter system. 

VI. CONCLUSIONS 

We have presented a model allowing for a systematic 
investigation of fusion barriers between heavy nuclei. It 
has been demonstrated that our approach which is based 
on a self-consistent semi-classical description of projec- 
tile and target nucleus and has no adjustable parameter 



predicts fusion barriers that are quite close to the ones 
obtained with the phenomenological proximity approach 
of Myers and Swiatecki ^7$ . We have given, in addition, 
an analytical form allowing for a simple approximation 
of the ETF fusion barriers, thus providing a very simple 
evaluation of fusion potentials for reactions throughout 
the periodic table. We believe that such an approach 
could serve as a guide-line for the research of the synthe- 
sis of super-heavy elements. 

It is also quite remarkable that our simple phenomeno- 
logical ansatz, Eq. nitty , adjusted to the reactions leading 
to super-heavy elements was able to reproduce even the 
fusion barrier of very light systems in a quite satisfactory 
way with a relative error of 2 % at most. 

One might also speculate about the importance of de- 
formation on the barrier heights obtained in our ap- 
proach. It seems immediately evident that taking into 
account deformation as an additional degree of freedom 
of the system target-projectile could often lead to a de- 
crease of the fusion barrier [23] . One knows, indeed, that 
e.g. the Coulomb barrier is substantially lower for two 
nuclei with ellipsoidal deformation and tip-to-tip orien- 
tation as compared to a different orientation or in the 
absence of deformation. It is obvious that for a system 
like 110 Cd + 154 Sm, which we briefly discussed above, 
where at least the target nucleus shows a substantial de- 
formation, such an effect plays a non negligible role which 
implies that taking deformation into account in our ETF 
approach might lead to much lower fusion barriers for 
those systems. 
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